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Acoustic Oscillations in Solar and Stellar Flaring Loops 

V.M. Nakariakov 1 , D. Tsiklauri 1 '*, A. Kelly 1 , T.D. Arber 1 , and M.J. Aschwanden 2 



^1" 

o 
o 

(N 



> 

(N 
(N 
(N 
O 

o 

Or 

6 ■ 
a: 



1 Physics Department, University of Warwick, Coventry, CV4 7AL, UK 

2 Lockheed Martin, Advanced Technology Center Solar & Astrophysics Lab, Dept. L9-41, Bldg. 252, 3251 Hanover 



Street, Palo Alto, CA 94304, USA 



Received ??? 2003 



Abstract. Evolution of a coronal loop in response to an impulsive energy release is numerically modelled. It is shown 
that the loop density evolution curves exhibit quasi-periodic perturbations with the periods given approximately 
by the ratio of the loop length to the average sound speed, associated with the second standing harmonics of an 
acoustic wave. The density perturbations have a maximum near the loop apex. The corresponding field-aligned 
flows have a node near the apex. We suggest that the quasi-periodic pulsations with periods in the range 10-300 s, 
frequently observed in flaring coronal loops in the radio, visible light and X-ray bands, may be produced by the 
second standing harmonic of the acoustic mode. 
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1. Introduction 

Wave activity of the solar corona attracts attention in 
relation with coronal heating and solar wind accelera- 
tion problems, and as an efficient tool for MHD coronal 
seismology (e.g. Nakariakov 2003). The observational ev- 
idence of coronal waves and oscillations is abundant. Low 
period coronal oscillations, in the range between a few 
seconds to several minutes, are believed to be associated 
with magnetohydrodynamic waves. In particular, propa- 
gating slow magnetoacoustic waves have been identified 
in polar plumes (Ofman, Nakariakov & DeForest 1999) 
and over loop footpoints (Nakariakov et al. 2000, Tsiklauri 
& Nakariakov 2001), and standing global slow modes in 
loops (Ofman & Wang 2002). The majority of confidently 
interpreted examples of the coronal wave activity has been 
found in the EUV coronal emission. 

Observations in other spectral windows, in particular 
in the radio band, also demonstrate various kinds of os- 
cillations (e.g., the quasi-periodic pulsations, or QPP, see 
Aschwanden 1987 for a review), usually with periods from 
a few seconds to several tens of seconds. It is commonly 
accepted that the waves with periods about several sec- 
onds may be produced by either propagating or standing 
sausage modes (Roberts, Edwin & Benz 1984; Nakariakov, 
Melnikov & Reznikova 2003, Aschwanden, Nakariakov & 
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Melnikov 2004). However, it is known that the existence 
of the non-leaky sausage mode cut-off imposes an upper 
limit on the possible wave periods: waves with the periods 
longer than the cut-off period cannot be trapped in the 
loops. For a magnetic cylinder of the radius a, the cut-off 
period P c is estimated (see, e.g. Roberts et al. 1984) as 
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where C s q is the sound speed in the loop, and C'aq and 
CA e are the Alfven speeds inside and outside the slab, 
respectively. For example, for typical flaring loop param- 
eters a — 6 Mm, Cao = 1 Mm/s, Ca = 3Cao and 
C s q = 0.5 Mm/s, the cut-off period is about 15 s, and 
it is difficult to make it greater than 20 s. However, X- 
ray band observations often give much longer periodici- 
ties, frequently in association with a flare. For example, 
periodicities from 20 s to 25 min have been presented by 
Harrison (1987), McKenzie & Mullan (1997) and Terekhov 
et al. (2002). Similar periodicities have been observed in 
the decimeter and microwave bands. In particular, Wang 
& Xie (2000) observed QPP with the periods of about 
50 s at 1.42 and 2 GHz (in association with an M4.4 X- 
ray flare). The coincidence of QPP periods observed in the 
X-ray and in radio bands is not a surprise, as the higher 
frequency radio bursts are found to correlate very well 
with X-ray bursts (e.g., Benz & Kane 1986). Moreover, 
pulsations with the periods significantly greater than the 
estimated sausage mode cut-off period have been found 
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in both hard X-ray and microwave bands simultaneously 
(e.g. Fu, Liu & Li 1996, Tian, Gao & Fu 1999). 

Similar oscillations, with the period of 220 s, have 
recently been observed in the white-light emission as- 
sociated with stellar flaring loops (Mathioudakis ct al. 
2003). Possibly, the 160 s periodic oscillations observed by 
Houdebine et al. (1993) during a flare on Ad-Leonis and 
26 and 13 s coherent oscillations observed by Zhilyaev et 
al. (2000) in EV Lac flares have the same nature. 

A possible interpretation of these medium period 
QPPs may be connected with kink or torsional modes 
(Zaitsev & Stepanov 1989). However, these modes are 
practically incompressible and, in the case of a small am- 
plitude, the produced perturbation of the magnetic field 
is also very weak. (E.g., the direct observations of kink 
modes in EUV, Aschwanden et al. 1999, Nakariakov et al. 
1999). Thus, it is not simple to link these weak perturba- 
tions of the magnetic field with observed QPPs. 

In this study, we suggest an alternative mechanism for 
the generation of long-period QPPs. We demonstrate that 
in a coronal loop an impulsive energy release generates 
the second spatial harmonics of the acoustic mode. These 
oscillations are of high quality and do not experience dis- 
sipation. This mode, producing noticeable perturbations 
of the loop density and generating field-aligned flows, is 
shown to be responsible for QPP with medium and long 
periods. 

2. The model 

We describe plasma dynamics in a coronal loop by solving 
numerically ID radiative hydrodynamic equations (infi- 
nite magnetic field approximation) that resemble closely 
the Naval Research Laboratory Solar Flux Tube Model 
(SOLFTM, Mariska 1987). The numerical code that we 
use is a ID version of Lagrangian Re-map code (Arber et 
al. 2001) with the radiative loss limiters. As in SOLFTM 
the coronal loop is connected with a dense, cold (10 K), 5 
Mm thick plasma region, that mimics the chromosphere, 
and which because of its large density provides sufficient 
amount of matter to fill the loop during the flare. The 
model includes: the effects of gravitational stratification, 
heat conduction, radiative losses, added external heat in- 
put, the presence of Helium, non-linearity, and Braginskii 
bulk viscosity. For the radiative loss function we use the 
form given in (Rosner et al. 1978) extended to a wider 
temperature range (Peres et al. 1982, Priest 1982). 

We performed a series of numerical runs simulating 
the response of the loop to a flare-like impulsive heat de- 
position at a chosen location. The heat deposition was 
modelled by the function 
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Fig. 1. Time evolution of the average temperature (upper 
panel) and the number density at apex (lower panel) in 
the flaring loop in the case of apex heating. The upper 
curves correspond to longer flares and the lower to the 
shorter flares. The vertical straight lines show the time 
interval of particularly distinct oscillations on the upper 



where, E$(s) and E^{t) were spatial and temporal parts 
of the heating function. Eq is the amplitude in units of 



erg cm -3 s ; sq is the position of the maximum heat 
deposition; a B is the heat deposition localization scale. 
Function J5J allows us to model spatial and temporal pro- 
files of the heat deposition separately. As it is discussed 
below, the results are weakly dependent upon the partic- 
ular choice of this function. The only restriction is the 
spatial and temporal characteristic scales to be shorter 
than the loop length and the acoustic transit time along 
the loop, respectively. 

In the numerical runs presented here, the size of the 
heat deposition region is a s — 7 Mm. The flare peak time 
is fixed in at 2200 s. The heat deposition duration, cr t , is 
variable, from 42 s to 333 s. The heating amplitude Q p is 
fixed at 2 x 10 4 , while E is fixed at 0.004 erg cm~ 3 s -1 . 
The loop length is taken to be 55 Mm. We keep the coronal 
part of the loop initially at a temperature of 1 MK and at 
a mass density of p = pm p n c = 6.6 x 10 -16 g cm -3 (at the 
loop apex), for a helium-to- hydrogen number density ratio 
of 0.05, i.e., with the mean molecular weight /i = 1.1. This 
corresponds to an electron density of n c = 3.6 x 10 8 cm -3 . 
In our approach, specific initial conditions in region con- 
necting corona to the chromosphere and chromosphere it- 
self have little or no effect on the corona dynamics as they 
are rapidly modified, self-consistently, according to the ra- 
diative hydrodynamics equations. Resolution in all numer- 
ical runs was fixed at 1000 grid points, which were dis- 
tributed non-uniformly in order to properly resolve strong 
gradients in the transition region. Even tripling the num- 
ber of grid points does not alter the numerical results, thus 
assuring full convergence of the simulation. 
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3. Acoustic second standing harmonics 

The typical loop evolution is shown in Figure ^ The im- 
pulsively heated loop responds with an increase of the 
temperature. The particular response is determined by the 
duration of the heat deposition er t - In Fig. ^ the cases of 
<r t = 42 s and 333 s are presented. The chromospheric 
evaporation increases loop density with up-flows of the 
order of a hundred km s _1 . Then, during the peak of 
the flare, the combined action of the heat input and the 
conductive and radiate losses, yields an oscillatory flow 
pattern with typical amplitudes of up to few tens of km 
s^ 1 . In the cooling phase, the matter drains out of the 
loop exhibiting oscillations; and the loop finally reaches 
an equilibrium. 

In the context of this study, the most interesting fea- 
ture in Fig. 1 is the clear quasi-periodic oscillations of 
the density time-profiles, especially about and after the 
peak of the flare. This behaviour closely resembles the 
radio and X-ray flaring light-curves (see the citations in 
Introduction). 

Typical time-distance plots of the evolution of the den- 
sity and the velocity are shown in Figure 2. Here, the slow, 
in comparison with the oscillation period, aperiodic vari- 
ations of the density and the velocity are subtracted. The 
velocity perturbation clearly demonstrates the node near 
the apex, while the density perturbation is maximal there 
(anti-node). This behaviour is consistent with the second 
standing acoustic mode, 
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where C s is a speed of sound, A is wave amplitude, L is 
loop length, and s is a distance along the loop with the 
zero at the loop top. 

The top panel of Figure |3| shows characteristic time 
signatures of the density and velocity oscillations at 
two points symmetric with respect to the loop apex. 
Depending upon the leg of the loop, the phase shift be- 
tween the density and the velocity oscillations is plus or 
minus quarter period, as it is expected for the second 
standing acoustic harmonic (see Eqs. J3J)). Perturbations 
of the density arc in phase at both points of measurement. 
The oscillations are quasi-harmonic, with a well-defined 
period. In the discussed case, the period of the oscillation 
is 64 s. According to Eqs. J2J), the period of the second spa- 
tial harmonic should be P — L/C s . The practical formula 
for the determination of the oscillation period is 



P/s « 6.7 x (L/Mm)/V(T/MK) 



(5) 



where T is the average temperature in the loop. (This 
expression is a particular case of the general formula sug- 
gested by Roberts et al. (1984), who pointed out that the 
second standing harmonic could likely be excited in coro- 
nal loops.) Substituting the loop length L = 55 Mm (see 
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Fig. 2. The time-distance plots of perturbations of the ve- 
locity and density, after subtraction of the slow evolution 
of the background, for the time interval highlighted for 
the upper curve in Fig 1. 




/BOO 



Fig. 3. Time signatures of the velocity and density pertur- 
bations outputted at the points —6 Mm (top panel) and 
+6 Mm (bottom panel) of the loop length for the same 
time interval as Fig El The solid curve shows the plasma 
number density in units of 10 11 cm -3 . The dashed line 



gives the velocity normalized to 400 km 



Here, the 



slow variations of the density and the velocity are sub- 
tracted. 



Fig. 2) and the average temperature, 25 MK (see top left 
panel in Fig. 1 in the range of 2500-2800 s), we obtain 
73.7 s, which is reasonably close to the result of the simu- 
lations. The discrepancy of about 15% is explained by the 
fact that the actual reflection points are situated over the 
loop footpoints (see Fig. |2J because of stratification. This 
reduces slightly the effective length of the loop and, conse- 
quently the resonant period of the mode. However, as the 
uncertainty in the observational determination of the loop 
length is comparable with this error, formula J5j provides 
a satisfactory estimation of the oscillation period. 
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4. Conclusions 

We suggest that the second standing acoustic mode may 
be responsible for QPP with periods of about 10-300 s (es- 
timated by Eq. (J3J) observed in solar and stellar flare light 
curves. This mechanism is similar to the interpretation of 
coronal loop oscillations observed with SUMER, proposed 
by Ofman & Wang (2002). The main new element in our 
study is connected with the mode excitation. We demon- 
strate that the second standing acoustic harmonic appears 
as a natural response of the loop to an impulsive energy 
deposition. The SUMER oscillations are likely to be as- 
sociated with some other excitation mechanism, as only 
a small fraction of SUMER oscillations are observed in 
association with solar flares (Wang et al. 2003). 

Traditionally, the acoustic wave interpretation was ex- 
cluded as these waves were supposed to be highly dissipa- 
tive. However, numerical simulations discussed above, as 
well as recently gained abundant observational evidence of 
the presence of acoustic waves in the solar corona, suggest 
that the observed periodicities can be associated with this 
mode. The physical mechanism responsible for the induc- 
tion of the quasi-periodic pulsations can be understood in 
terms of auto-oscillations generated by an electric-circuit 
generator. Indeed, the physical system modelled here con- 
tains all the necessary ingredients of a generator: the DC 
power supply (thermal instability), the nonlinear clement 
(the plasma) and the resonator (the loop). This may ex- 
plain why the oscillations may be observed to be dissipa- 
tionless. However, proper analytical theory of the excita- 
tion of this mode is still to be developed. 

The oscillation period is determined by the ratio of the 
loop length and the average sound speed in the loop. The 
typical amplitude of the density and velocity perturba- 
tions is 2-10% of the background. This value is consistent 
with observed amplitudes of X-ray QPPs (McKenzie & 
Mullan 1997). However, the observed amplitudes of ra- 
dioband QPPs are sometimes higher. This discrepancy 
would be resolved by taking into account the specific 
mechanism responsible for the modulation of radio emis- 
sion, and the line-of-sight angle (e.g. Cooper, Nakariakov 
& Williams 2003). In certain condition, the modulation 
mechanism can amplify the pulsations up to the observ- 
able level. 

The model developed here is quite approximate, as it 
does not take into account two-dimensional MHD effects 
such as centrifugal force and the perturbation of the loop 
cross-section. However, we believe that the neglected ef- 
fects do not change the qualitative picture described here 
and that Eq. (JSJ) gives a correct estimation for the oscilla- 
tion period. 

Finally, we would like to emphasize that the genera- 
tion of the second spatial harmonic of an acoustic wave, 
by means of which we try to explain the observed quasi- 
periodic oscillations in flaring loops, is a consistent feature, 
seen for a wide range of physical parameters, including the 
case of non-symmetric heating functions. More detailed 
analysis will be presented elsewhere. 
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